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Boundary losses and spatial statistics of complex modes in a chaotic microwave cavity 
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We experimentally study the various manifestations of ohmic losses in a two-dimensional mi- 
crowave chaotic cavity and exhibit two different contributions to the resonance widths. We show 
that the parts of these widths, which vary from mode to mode, are associated to ohmic losses located 
at the boundary of the cavity. We also describe how this non-proportional damping is responsible 
for the complex character of wavefunctions (corresponding to a spatially non- uniform phase) , which 
is ubiquitous in open or dissipative wave systems. We experimentally demonstrate that the non- 
proportional widths are related to a single parameter, which measures the amount of complexity of 
wavefunctions, and provide theoretical arguments in favor of this relation. 
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Investigations of the morphology of wavefunction in 
partially open chaotic wave systems have recently be- 
come an active domain of research, especially through 
theoretical works concerned with quantum transport in 
the ballistic regime0,|3|- That the spatial properties of 
resonances are closely related to their spectral widths 
has been recognized ever since the pioneering work of 
Porter and Thomas Q until more recent experimental in- 
vestigations in a superconducting microwave cavity 
These experiments resort to a chi-squared distribution 
of widths based on results of Random Matrix Theory 
(RMT), which is known to describe spectral and spatial 
statistics of closed chaotic wave systems 0. Microwave 
cavities turn out to be ideal model systems to study 
open wave systems 0, since they offer the possibil- 
ity of opening channels through antennas or connected 
waveguides. They can also be qualified as open be- 
cause of dissipation through ohmic losses at their bound- 
aries. Locally, boundaries may be viewed as distributed 
equivalent opening channels with ad hoc impedances. 
Through such losses, the standing waves of a closed cav- 
ity are changed into a mixture of standing and traveling 
waves, the latter being associated to the energy flow- 
ing out of the systemjsj. Proportional damping, i.e. 
the same decay rate for all modes in a given frequency 
band, can only lead to modes that are real (up to a 
global phase) and associated to the standing wave part 
of the wavefunction. Hence the observation of complex- 
valued wavefunctions implies non-proportional damping 
mechanisms giving birth to the traveling part of the 
wavefunction|22j . A thorough description of the link be- 
tween non-proportional damping and the complex char- 
acter of wavefunctions is still lacking, as recently empha- 
sized by Lobkis and Weaver in the context of acoustical 
reverberant dissipative system^. 

In this letter, we present experiments in two-dimen- 
sional (2D) chaotic microwave cavities at room tempera- 
ture where dissipation is essentially due to ohmic losses 



at boundaries. In such chaotic cavities, modes are gener- 
ically ergodic, i.e. present speckle-like intensity pat- 
terns. The importance of the imaginary part of the 
modes implied by the presence of losses can be measured 
through a single statistical parameter, namely the ratio 
{{Imip)'^) / {(Reip)'^) = where the wavefunction ip is 
conventionally assumed to be real {q = 0) in the case of 
vanishing losses. The principal aim of this paper is to 
demonstrate the intimate relation between this parame- 
ter q and the part of the decay rate (or modal width) asso- 
ciated to non-proportional damping. We provide strong 
experimental evidences of this relationship by carefully 
analyzing hundreds of resonances in a Sinai-like chaotic 
cavity. 

The cavity is the slice of a rectangular waveguide closed 
at both ends, with contour C of length L = 2.446 m, sec- 
tion S of area A = 0.352762 m^ and thickness d = 5 mm. 
Inside this flat parallelepipedic cavity is inserted an ob- 
stacle with the shape of a disk of radius 6 cm and thick- 
ness d. As long as the wavelength A = 27rfc~-'^ is larger 
than d, the cavity may be considered as two-dimensional 
(2D) and only admits transverse magnetic modes. If the 
conductor is perfect, the electric field E = E^'^\x,y)z is 
perpendicular to the plane (x, y) of the cavity, the mag- 
netic field is parallel and there is no dissipated power. 
This ideal situation will correspond to the zeroth or- 
der of our description of the field near the surface of 
the actual conductor. Denoting ipQ = E'^^'>{x,y), the 
time- independent Maxwell equations are reduced to a 2D 
Helmholtz equation: 



(V' + e/xw2)V'o = , 



(1) 



where e and fi are respectively the permittivity and the 
permeability inside the cavity. On the contour C, ipo 
obeys Dirichlet conditions ipQ = 0. 

Ohmic losses at the contour are easily computed by 
following a boundary perturbation technique jlfll| . Indeed, 
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the perturbed boundary conditions (BC) now read: 
^~-(l + i)^a„V onC, 



(2) 



where 5 = ^J2/ ^c<^contUJ is the skin depth (/ic being the 
magnetic permeabihty of the conductor and cJcont the 
effective conductivity of the contour), and where 9„ de- 
notes n • V (n being the unit normal vector directed to- 
ward the interior of the conductor). Green's theorem 
apphed to this 2D problem straightforwardly yields: 



(3) 



where oj is the perturbed eigenfrequency, ojo being the 
corresponding real eigenfrequency of the lossless problem. 
Writing u) = loq + 6lo — iV equation Q becomes: 



(5cj-ir/2~-(l+i) 



1 yUc(5 /^d£|i9„Vl 



(4) 



leading, in the case at hand, to equal corrections on both 
real and imaginary parts of lo. In our context fJ-c/lJ- is 
practically equal to unity. A similar perturbative ap- 
proach can be applied to the ohmic losses at the top and 
the bottom of the cavity. One finds that they lead to ho- 
mogeneous proportional damping. The total width due 
to ohmic losses can thus be split into two parts[ll|: 



pf2 _ pP 



L 



where ^ is defined by: 



(5) 



(6) 



The first term in equation © is the proportional 
(homogeneous) damping part of the width (depending 
only on uj without any dependence on the spatial struc- 
ture of the mode) due to effective conductivity aends of 
the top and bottom. The second term is the non- 
proportional part of the width related to a parameter ^, 
which depends on the spatial structure of the mode at 
the contour, and therefore varies from mode to mode. 

Here it is particularly interesting to revisit equation (0) 
which turns out to give a quite natural hint of how ohmic 
losses on the contour induce a small amount of imaginary 
part for wavefunctions that are purely real in the unper- 
turbed limit. Indeed, assuming an incident plane wave 
with unit amplitude, the BC (0) lead to a reflected wave 
with a small (of order kS) dephasing. The same BC lead 
to equation which relates the width to the quantity 

For the rectangular cavity the values of ^ are read- 
ily calculated and display fluctuations around unity, as 



generally expected[10j. Thus we could check the validity 
of formula 0) in our experimental situation|l2j. In the 
case of a Sinai-like chaotic cavity one may view a typical 
speckle-like mode locally as the superposition of plane 
waves with random directions, amplitudes and phases. 
This led Berry[l3 to conjecture that the wavefunctions 
should be ergodic on the average, i.e. should display all 
statistical features of a Gaussian random field. Apart 
from specific non-ergodic features like scars (local inten- 
sity enhancement along short unstable periodic orbits), 
the relative fluctuations of ^ about unity are expected to 
decrease with the wavelength A as the contour C encom- 
passes a larger number ~ L/(A/2) of equivalent opening 
channels. 

The cavity used for the measurements presented here is 
made of OFHC Copper. Ten antennas are fixed through 
the bottom plate and are used as emitters and receivers 
as well. These antennas (diameter 1.27mm) are only 
2 mm inside the cavity. A rather low coupling is thus 
achieved, making the role of the antennas nearly negli- 
gible regarding their scattering and loss characteristics. 
Radiation losses are also made negligible by a tight fitting 
of the various parts of the cavity. Losses are therefore 
dominated by ohmic losses on the walls of the cavity. A 
Vector Network Analyzer HP 8720 D is used for the mea- 
surements. We restrict our studies to frequencies rang- 
ing up to 5.5 GHz since, at higher frequencies, the modal 
overlap defined as F x n{uj) {n{uj) being the modal den- 
sity per unit interval of lo) becomes too important and 
invalidates any perturbative approach. Below 5.5 GHz, 
apart from few broader resonances, the modal overlap 
remains smaller than unity. The measurements are per- 
formed by 0.5 GHz-bands, sampled with a frequency step 
0.3125 MHz. The 45 possible couples of antennas are 
used. On the whole, it amounts to 756000 points mea- 
sured between 0.25 GHz and 5.5 GHz. 

In order to fit the transmission between two dis- 
tinct antennas, we use a Breit-Wigner formula (see e.g. 
reference 0) which, in our case, reads [13 



i^riFj^ 



The Lorentzian associated to the nth resonance for 
the transmission between the antennas a and b is thus 
characterized by its complex amplitude a'^ -f i/3^'' = 
~2iT'^{LUn)i^n['ra)tpn{'r'b), its width F^^ and its central 
frequency uon = ~ Tn/^j where a;° stands for the 
jith eigenfrequency of the unperturbed cavity. T{lj) 
is a function representing the coupling of the antennas 
(all supposed to be identical) to the cavity whose ex- 
plicit form, which essentially depends on the geometry of 
the antenna, has been derived in^^l- We were able to 
test our model for antenna coupling very precisely in an 
empty rectangular cavity where wavefunctions and eigen- 
frequencies are readily calculated. For the present pur- 
pose, it is enough to know that T is a slowly varying 
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function on a frequency scale much larger than the mean 
spacing n~^{uj). This coupling is obviously also respon- 
sible for additional losses and can be accounted for by 
introducing corresponding partial widths. The contribu- 
tion of the latter to the total widths being at most of the 
order of one percent for all the measurements presented 
here, we could safely neglect them. 

Assume that the real and imaginary parts of the wave- 
function '0(r) are Gaussian centered random variables 
with standard deviations a,. = ((Re ip)'^) ^'iid (Ti = 
^ ((Im V')^)- The key parameter q = ai/ar is a mea- 
sure of the openness of the system. Its value is when 
the cavity is completely closed (real wavcfunctions) and 
tends to unity for totally absorbing boundaries (complex- 
valued wavcfunctions with equal ar and ct^). In practice, 
as the spatial sampling over ten antennas for a given res- 
onance is undersized to estimate ar and ai accurately, 
we are led to restrict ourselves to a coarse-grained esti- 
mation of q by grouping several contiguous resonances 
in frequency windows. This may be achieved by build- 
ing the probability distribution of the phase ip"^^ of the 
complex amplitude of a"'' + 1/3°'': 



(8) 



where tpn (r) is the phase of the complex field {"t) ■ This 
distribution is readil y o btained from the Gaussian hy- 
pothesis and reads |9l ll-d llq: 
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2 + (l-g2)2sin^(^a6^|) 



(9) 



This distribution displays a peak centered on = 7r/2 
(see figure The width of the peak increases with 
q so that the distribution tends to a Dirac delta as 
q ^ or to a uniform distribution as g 1. The 
distributions (a) and (b) shown in figure ^ have been 
obtained by using 30 resonances lying in the frequency 
intervals [3.76 GHz, 4.1 GHz] and [5.24 GHz, 5.5 GHz] re- 
spectively, yielding coarse-grained estimations of q equal 
to 5.0 X 10~2 and 10.1 x lO"^ respectively. 

In order to establish the connection between the spatial 
statistics of complex wavcfunctions and their (spectral) 
widths due to non-proportional damping, it is useful to 
introduce a ray picture of losses due to reflections at the 
contour. For a plane wave incident on an air- metal inter- 
face with angle of incidence 6*, it can be shown that 
the amount of absorbed energy reads: 



4a cos 6 



2a2 cos2 + 2acos6' + l ' 



(10) 



where a ~ ujej {^acxmt)- In the case of OFHC copper 
a ~ 10"^ around 1 GHz and the previous equation re- 
duces to 




FIG. 1: Experimental histograms of phase yp" ob- 
tained with 30 resonances in the frequency intervals (a) 
[3.76 GHz, 4.1 GHz] and (b) [5.24 GHz, 5.5 GHz]. Fitted dis- 
tributions @ in solid line: (a) q = 5.0 x 10~^ and (b) 
q = 10.1 X 10"^ 



For an ergodic wavefunction in a chaotic billiard, the 
appropriate statistical measure for angular average is 
d(sin Q) /2 leading to a mean absorption coefhcient 



{T) = 'n:\Juje/{2acont) ■ 



(12) 



Having in mind the analogy of our problem with ultra- 
sonics in a reverberant body^3l or room acoustics [l^. 
one may use the so-called Sabine's law[l8j to deduce the 
average width Trefi related to absorption at the contour 



r.e/;-^^(r)-^y^|^ 



(13) 



This last result may be identified with the expression 
of F^-^ in equation jSJ provided that ^ be replaced by 
its mean value for an ergodic mode namely unity [ll|. 
As q and the transmission T both vary within the same 
interval ranging from (completely closed system) to 1 
(completely transparent contour) , wc propose to identify 
q and T, leading to the following relationship: 



cL 

— 7( 



(14) 



T ~ 4a cos . 



(11) 



This simple identification is all the more natural that the 
spectral perturbation F^'^ and the spatial perturbation 
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q are of the same order. Indeed, the relationship ^ 
indicates that q is of order k5. Likewise, equation (0) 
shows that {A/cL) Y'^^ is also of order k5. This excludes 
any power law relationship, other than proportionality, 
between q and T^^ . 

From our measurements, we used formulas Q and 
^ to extract, by an appropriate and original fitting 
procedure^3|j the values of the non-proportional widths 
pATP Jqj, gT^Y regonances up to 5.5 GHz. These data 
are shown on figure [21 where an average over contigu- 
ous groups of 15 (a) and 10 (b) resonances is performed. 
They are compared to expression (|14|l where a coarse- 
grained evaluation of q is deduced from the distribution 
© built with the same groups of resonances. With only 
10 resonances, we have checked that fitting the experi- 
mental histograms of phase with formula jnj is still le- 
gitimate. The obvious correlation between the effective 
measure of imaginary parts, given by the parameter g, 
and the variations of the widths gives strong support to 
our hypothesis embodied in relation H14|l. The agreement 
is fairly good even though some small discrepancies can 
be observed, which can be attributed to the lack of er- 
godicity of certain modes. Our behef is that the relation 
((TUl should hold for each truly ergodic mode but its ver- 
ification is out of reach with our present experimental 
set. 
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FIG. 2: Solid lines correspond to evaluations of extracted 
from our measurements through expression Q with an aver- 
age over contiguous groups of: (a) 15 resonances and (b) 10 
resonances. Dashed lines shows relation 11411 with coarse- 
grained q deduced from the distribution 0. 



That the wavefunctions of a lossy wave system are in 
general complex has not been well appreciated in the past 
literature maybe because the most common and simple 
model of losses leaves the modes real. Such a model pre- 
dicts that the resonance widths should only vary slowly 
with frequency. In this letter, for the first time to our 
knowledge, we propose theoretical arguments in favor of 
the hypothesis that the imaginary parts of the modes 
should correlate with the non-proportional part of the 
widths. Through an original analysis of transmission 
measurements in a lossy chaotic microwave cavity, we 
provide clear experimental evidence of its pertinency. It 
should be stressed that our point of view is complemen- 
tary of RMT in that it focuses on spatial aspects whereas 
RMT for open systems is mainly concerned with spectral 
statistics|5|- We believe that the present work may have 
conceptual implications in transport in open mesoscopic 
devices Il9l and also in structural acoustics l23- 
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